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ABSTRACT 

The optical lightcurve of GRB 010222 exhibited one of the slowest decays of any gamma-ray burst 
to date. Its broadband properties have been difficult to explain with conventional afterglow models, as 
they either require the power law index of the underlying electron energy distribution to be low, p < 2, 
or that the outflow is quasi-spherical thus reviving the energy problem. We argue that the slow decay 
of GRB 010222 and a linear polarization of 1.36 ± 0.64%, is naturally explained by a jet model with 
continuous energy injection. The electron energy distribution then has p = 2.49 ± 0.05, fully consistent 
with the expectation from detailed modelling of acceleration in relativistic shocks, that p > 2, thus 
alleviating the "p-problem" . 

Subject headings: gamma rays: bursts — polarization 



1. INTRODUCTION 

GRB 010222 was a bright gamma-ray burst (GRB) lo- 
calized by BeppoSAX 

(Piro et al. 2001). X-ray observations were reported by 
in 't Zand et al. (2001). The optical afterglow was dis- 
covered by Henden & Vrba (2001), 4.3 hours after the 
burst and a redshift of z = 1.477 was determined by 
Jha et al. (2001). Further optical/near- infrared obser- 
vations have been reported by Stanek et al. (2001), Lee 
et al. (2001), Masetti et al. (2001), Cowsik et al. (2001), 
Sagar et al. (2001) and Mirabal et al. (2002). In our anal- 
ysis we shall adopt the light curve fit of Henden et al. 
(2002), that is based on data extending up to 80 days 
after the burst. Their best fit in the i?-band gives a pre- 
break power-law slope of ai = —0.66 ± 0.03 and a post- 
break slope of a2 — —1.40 ± 0.02, with a break time of 
tb = 0.58 ± 0.04 days. The observed optical spectral in- 
dex is about P = —1.0 (e.g. Stanek et al. 2001), but as it 
may be strongly affected by host extinction, it can lead to 
ambiguous inference of the afterglow properties. In par- 
ticular, the relationship between the light curve decay in- 
dices, the a's, and the spectral index, /3, that is predicted 
by synchrotron models of afterglows (e.g. Sari, Piran & 
Narayan 1998; Sari, Piran & Halpern 1999; Price et al. 
2002), can lead to inconsistent interpretations due to the 
unknown host extinction. 

To date, interpretation of the GRB 010222 afterglow 
observations has mainly relied on two possible scenar- 
ios: i) A narrow sideways expanding jet propagating in 
a low-density medium (e.g. Stanek et al. 2001). ii) A 
wide jet or spherical fireball transiting from relativistic 
to non-relativistic regime (e.g. Masetti et al. 2001; in 't 
Zand et al. 2001). The latter models require a very dense 
medium and implies a very large energy release, thus 
bringing back the energy problem (Kulkarni et al. 1999; 
Andersen et al. 1999). Both classes of models require a 



hard electron energy distribution, i.e. p < 2, in the first 
case even as low as p ~ 1.4. This is contrary to most 
other bursts that seem to be adequately fit with models 
where p « 2.3 — 2.5 (van Paradijs, Kouvelioutu & Wijers 

2000) . Recently, Panaitescu & Kumar (2002) presented 
fits to several afterglows, concluding that p is smaller 
than 2 in a number of them. Such small inferred values 
of p signal a departure from the standard fireball model 
(Meszaros 2002), and introduce additional free parameters 
into the model, such as an upper cutoff in the electron 
energy distribution, or additional assumptions about the 
electron acceleration mechanism to facilitate generation of 
flat energy distributions (Dai & Cheng 2001; Bhattacharya 

2001) . The p-problem then arises from the fact that 
detailed modelling of particle acceleration in relativistic 
shocks indicates that p « 2.2 — 2.3 (Achterbeg et al. 2001). 

In addition, if p is assumed to be constant through- 
out the fireball evolution, the magnitude of the observed 
light curve break in the case of GRB 010222, Aa = 
oi\ — cti = 0.74 ± 0.04, cannot be explained by the above 
models. The sideways expanding jet model predicts a 
break magnitude of Aa = 1 — ai/3 = 1.22 ± 0.01, while 
the fireball transiting to the non-relativistic regime gives 
Aa = (ai -I- 3/5) = 0.06 ± 0.03 for slow cooling elec- 
trons. For fast cooling electrons the latter model predicts 
Aa = (ai + 1) = 0.34 ± 0.03. The break magnitudes pre- 
dicted by the models are in all cases very different from 
the observed value that is however, in perfect agreement 
with the prediction of a jet model with a fixed opening an- 
gle (Meszaros & Rees 1999), Aa = 3/4, being essentially 
of geometrical origin. 

In this Letter we use polarization measurements to ar- 
gue against interpreatations based on spherical models 
transiting to the non-relativistic regime. Furthermore, we 
demonstrate that the observations of GRB 010222 can be 
naturally interpreted with a jet model with a small open- 
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ing angle and continuous energy injection. 

2. OBSERVATIONS 
2.1. Polarimetry 

Polarimetric imaging observations were obtained at the 
2.56-m Nordic Optical Telescope (NOT) on La Palma, The 
Canary Islands, with the Andalucia faint object spectro- 
graph (ALFOSC), using two calcite plates together with a 
y-band filter. Each calcite plate provides the simultaneous 
measurement of the ordinary and the extraordinary com- 
ponents of two orthogonally polarized beams. Thus, one 
image gives either the 0° and 90° components or the 45° 
and 135° components simultaneously. The calcite plates 
produce a vignetted field of about 140" in diameter. The 
detector was a thinned Loral 2048 x 2048 pixel Charge- 
Coupled Device (CCD) giving a pixel scale of 0.188". The 
ordinary and extraordinary beams appear as images sep- 
arated by about 15" on the CCD detector. 

Observations began 21.88 hours after the burst with ob- 
serving and reduction procedures similar to that used for 
GRB 990123 (Hjorth et al. 1999), each image providing 
fluxes in two of the four orientations required. Individual 
exposure times were 600 or 900 s resulting in total ex- 
posure times of 2400 s per orientation at a mean epoch of 
Feb 23.25169 2001 UT, 22.65 hours after the burst (see Ta- 
ble 1). No information on the position angle is available, 
as no polarization standards were observed. The images 
were flat fielded using standard procedures with appropri- 
ate dome flats in each of the two orientations. The three 
images in each orientation were combined. Absolute pho- 
tometric calibration was performed relative to star A of 
Henden et al. (2002) which has V = 14.582 ± 0.022. Mete- 
orological conditions were fine and the measured resolution 
of the images was ~ 0.8-1.0". Figure la shows an excerpt 
of the combined images. 

Aperture photometry was performed on the combined 
images with the DAOPHOT II/ALLSTAR (Stetson 1987, 
1994) software package. Point-spread function (PSF) pho- 
tometry turned out to be difficult due to the lack of suit- 
able PSF stars and the strongly orientation-dependent 
PSF shape and width. Aperture growth curves were com- 
puted for the OA, star B and a galaxy (G) in the field (see 
Fig. la). Star A was saturated in these images. Aper- 
ture photometry of star A was therefore performed on the 
short exposures obtained immediately after the six long 
exposures (cf. Table 1). 

In order to correct for possible interstellar polarization, 
polarization induced by the telescope and instrument, and 
the significant variations in stellar shape as a function of 
orientation, we computed aperture growth curves for OA, 
B and G relative to A. 

The Stokes vector, the instrumental Q and U values, 
and the linear polarization were computed from the de- 
rived fluxes in the four orientations 

(Hjorth et al. 1999; Wijers et al. 1999; Rol et al. 2000) as 
a function of aperture radius. The results are dependent 
on the aperture radii used. For small radii the results are 
dominated by photon noise and small variations in the 
shapes and widths of the objects relative to star A. For 
large radii the noise from and difficulty in subtracting the 
sky background dominates. Thus, one can only hope for 
reliable results for intermediate aperture radii. We used 



two criteria for selecting the aperture radius, namely (i) 
minimizing photon noise which favored aperture radii be- 
tween 4 and 6 pixels and (ii) insensitivity to the aperture 
radii of the Q and U values which favored aperture radii 
between 5 and 7 pixels. We note that B never reached 
independence of aperture radius and shows a variation in 
Q larger than that expected from photon noise. We there- 
fore consider the results for B unreliable, probably due to 
relatively high photon noise from the fainter star. Aper- 
ture radii between 5 and 7 pixels for the OA and G are 
shown as filled circles in Fig. lb. The mean of these give 
P{OA) = 1.50 ± 0.64% and P{G) = 0.4 ± 0.9%. We note 
that the polarization of the galaxy is consistent with zero 
as it should if star A is unpolarized. 

The effects of depolarization or interstellar polarization 
in the Milky Way are expected to be negligible at the 
high latitude of GRB 010222 {b = 60.9°) - well out of 
the Galactic plane. For the low Galactic extinction of 
E(B-V) = 0.023 (Schlegel, Finkbeiner & Davis 1998), the 
maximum interstellar polarization or depolarization is less 
than 0.2% (Serkowski, Mathewson & Ford 1975; Berdyu- 
gin & Teerikorpi 1997). Thus the interstellar polarization 
or depolarization towards GRB 010222 is negligible. 

At low significance levels of the computed afterglow po- 
larization, a correction must be applied to account for the 
non-Gaussian nature of the underlying probability distri- 
bution (Wardle & Kronberg 1974). When corrected for 
this effect the V-band linear polarization of the optical 
afterglow of GRB 010222 is 1.36 ± 0.64%, similar to the 
degree of polarization observed in a number of other bursts 
(e.g. Covino et al. 1999; Rol et al. 2000). 

As seen in Table 1, the optical afterglow decayed dur- 
ing the observations. As a check on the reliability of our 
polarization measurements we computed the V-band mag- 
nitudes by taking the mean of the two orientations at each 
epoch. The resulting decay slope of —1.19 ± 0.25 is con- 
sistent with the observed lightcurve (e.g. Henden et al. 
2002). 

2.2. X-rays 

We have also analyzed data from Chandra X-ray Ob- 
servatory observations starting about 15 hours after the 
burst, yielding a net exposure time of 29.5 ksec. Due to the 
brightness of the X-ray afterglow the data are severely af- 
fected by pile-up. Hence, the spectra were analyzed follow- 
ing the procedure suggested by Davis (2001), taking pile- 
up effects explicitly into account in the spectral modeling. 
We extracted a spectrum from a circular 4 arcsec diameter 
region centered on the X-ray afterglow and fitted the 0.5- 
10 keV spectrum with three model components: i) An in- 
trinsic power law, ii) Galactic absorption fixed at the nomi- 
nal value n_f/ = 1.6x10^° cm~^ (Dickey & Lockman 1990), 
and iii) intrinsic absorption in the GRB host fixed at the 
anticipated redshift of z = 1.477. This model provides 
an excellent fit (x^ = 116 for 204 d.o.f.), with a best 
fit spectral index /3 = —0.72 ±0.17 and intrinsic GRB 
host absorption rijy = 6.5 ± 0.11 x 10^^ cm~^ (Itr errors). 
There is no evidence for additional spectral features. The 
4 keV flux at the mean epoch of observations, 0.81'^, was 
0.268 ± 0.054 /LtJy. The near contemporaneous BeppoSAX 
observations give /? = -0.97±0.05 (in 't Zand et al. 2001). 
The difference between our spectral index and the Bep- 
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poSAX results is due to diflFerent derived values of nn- 
Using their method and fixing hh at their value, we find 
/3 = -1.13 ±0.20. 

3. A JET WITH CONTINUOUS ENERGY INJECTION 

Sources of synchrotron radiation arc generally expected 
to exhibit polarization in their emission, with the degree 
of polarization as high as 60 — 70% (e.g. Rybicki & Light- 
man 1985). This is, however, strongly dependent on the 
degree of regularity in the magnetic field. In fireball mod- 
els, the field is expected to be highly entangled with no 
preferred direction, hence little or no polarization is ex- 
pected. This would also be the case in a spherical fireball, 
due to symmetry, even if the magnetic field had a regular 
component. 

In a jet-like fireball, if the magnetic field has a regular 
component and the line of sight is at an angle with the 
center of coUimation, some polarization may be observed 
(Ghisellini & Lazzati 1999; Sari 1999). The degree of po- 
larization observed in GRB 010222 thus suggests that the 
fireball is slightly asymmetric. It follows that a transition 
to a non-relativistic regime of a spherical fireball or a wide 
jet, is not likely to be a correct description of this event. 

On the other hand, the observed light curve break mag- 
nitude, Aa = 0.74 ± 0.04, is exactly what is expected in 
a jet of fixed opening angle (Meszaros & Rees 1999). We 
will henceforth adopt that geometry. 

The isotropic energy release of GRB 010222, as esti- 
mated from the SAX data is E52 — 154.2 ± 17.0 (in units 
of lO'"^^ erg), in a cosmology with fim = 0.3, Ha = 0.7 and 
ho = 0.65 (Amati et al. 2002). The initial Lorentz- factor, 
estimated from equation (10) in Sari & Piran (1999), is 
then Tq « 350, ignoring the weak dependence on the am- 
bient number density, no . Using the light curve break time 
tb = 0.58±0.04 days (Henden et al. 2002), we find that the 
corresponding coUimation angle at tb is ^0 = ^/^b ~ 2.0°, 
reducing the estimate of the energy released in the burst 
to E Ki 10^^ ergs, fully consistent with the results of Frail 
et al. (2001), although near the higher end of their distri- 
bution. 

We advance the hypothesis that the exceptionally 
slow light curve decay of GRB 010222 is the conse- 
quence of continuous energy injection. If the fireball 
energy injection is continuous over an extended period 
of time rather than instantaneous, it may be modeled 
as a power law in ejected mass (Rees & Meszaros 1998; 
Sari & Meszaros 2000). The resulting time evolution 
of the Lorentz-factor can then be expressed as F = 
^oit/to)~^^~^^^^'^~^^\ where s is the power law exponent 
of the ejected mass distribution and g is the power law 
index of the ambient density distribution. We will hence- 
forth assume a homogeneous environment (g = 0). Then 
s = 1 reproduces the instantaneous case as the en- 
ergy release is then constant, while s > 1 implies that 
the energy is dominated by material with low Lorentz- 
factors (Rees & Meszaros 1998). We will take s ~ 2 
as a representative case of moderate energy injection 
(Sari & Meszaros 2000). Knowing the hydrodynamic evo- 
lution allows one to derive the corresponding light curve 
and spectral evolution as in the standard fireball model 
with instantaneous energy injection. General expressions, 
valid also for inhomogeneous external medium, can be 
found in Sari & Meszaros (2000). We will use their no- 



tation in what follows. It should be emphasized that a 
reverse shock is expected to produce substantial flux at 
low frequencies, especially at late times. 

Interpreting the optical and X-ray light curves as due 
to the forward shock we find, using ai = —0.66 ± 0.03 and 
s = 2, that 



P = 



3{s+l) - (7 + s)ai 



2.49 ±0.05, 



(1) 



in agreement with many other bursts. It follows that the 
intrinsic spectral index is /? = ~{p — l)/2 = —0.75 ± 0.03. 
This is fully consistent with the X-ray observations, if the 
cooling frequency, Vc was above the X-rays at the time 
of observations, but it requires moderate extinction in the 
optical (Lee et al. 2001; Masctti et al. 2001). It provides a 
consistent picture of the broadband spectrum from optical 
through the X-rays. The uncorrected F-fiux at O.Sl'* is 
45/xJ, using the temporal index a = —1.19 as inferred from 
our polarimetry. Adopting the SMC extinction correction 
of Lee et al. (2001), we find that Ay = 0.3 corresponds to 
P = —0.72, consistent with the intrinsic spectral slope we 
inferred above and our X-ray observations. 

We note that if 2 < p < 3, when s = 2, then the pre- 
break slope would be in the range 0.33 < —ai < 1.0. For a 
jet with fixed ^o, we would then have 1.08 < —02 5! 1-75. 
It is not straight forward to estimate the post break slope 
in an expanding jet model, as it requires a knowledge of 
how the continuous energy injection modifies the jet evolu- 
tion subsequent to the moment when T < 1/6*0. This is the 
subject of a separate paper (Bjornsson et al. in prepara- 
tion). We expect, however, a tranistion period where the 
temporal slope depends both on p and s. The asymp- 
totic decay slopc^ will in the end approach — p, but this 
may not be reached in all cases before the optical tran- 
sient fades below detectability or the host galaxy starts to 
dominate the emission. The late time decay slope derived 
by Fruchter et al. (2001), a2 = -1.7 ± 0.05, may indeed 
indicate that the afterglow continued to steepen at very 
late times. Wc suggest that the above scenario may also 
apply e.g. to GRB 020813 as its light curve decay is very 
similar to that of GRB 010222 (e.g. Bloom, Fox & Hunt 
2002; Gladders & Hall 2002). 

The flux from the forward shock has a maximum at the 
frequency v^^, that in this interpretation is below the op- 
tical. The spectral region from optical through X-rays is 
therefore represented by the same power-law, F^, oc z/~° ''^. 
From our X-ray data at 0.81*^, we obtain 



1/2 7/8 3/2 

'■0 ^B,-2^e,0.5 



2.53-10-^. 



(2) 



To solve for the model parameters, we need two more con- 
straints. The optical and the X-rays are degenerate in 
the sense, that they are both between and 1/^. The 
reverse shock is expected to produce substantial emission 
around the frequency f^, that is deflned analogously to 
f^, and is related to it by = vl^/V"^, the peak flux of 
the reverse shock bcnng a factor F larger than the peak 
flux of the forward shock. For ~ 10^^ Hz, could be 
from a few hundred GHz to 10^^ Hz. Furthermore the flux 
from the reverse shock may be expected to stay constant 
or even increase for a period of time (Panaitescu, Meszaros 
& Rees 1998), or decay as t'^/^ (Sari & Meszaros 2000). 
The radio and sub- mm observations of Frail et al. (2002), 
show a constant flux onwards from about 0.35*^ and 0.8*^ 
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at 350 and 250 GHz, respectively. They convincingly ar- 
gue that this flux originates in the host galaxy rather than 
the afterglow. That sets an upper limit on the emission 
predicted by the model, in principle enabling us to further 
constrain the parameters. In practice, as jy^ is unknown, 
the constraints in this case are very weak. It is apparent 
though from eq. (2), that for reasonable values of no, and 
ee,o.5) the magnetic field is very weak. 

A generic property of extended injection models is that 
the afterglow fades more slowly than in the case of in- 
stantaneous injection. We have chosen a simple power-law 
distribution of T as an example of such a scenario. Other, 
perhaps more detailed models, are of course possible (e.g. 
Zhang & Meszaros 2002), but would only affect the above 
interpretation in the details. The interpretation presented 



here, requires the introduction of an additional model pa- 
rameter (s), but this is also true in other models for the 
case of p < 2 as discussed above. An inferred value of p > 2 
is furthermore supported by detailed modelling of particle 
acceleration in relativistic shocks (Achterbeg et al. 2001). 
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Table 1 

Journal of polarimetry observations of GRB 010222 



Time orientation exp. time V 

(Feb 2001 UT) (deg) (s) (mag) 

23.22499 45/135 900 20.423 ±0.011 

23.23640 0/90 900 20.444 ±0.011 

23.24782 45/135 900 20.436 ±0.011 

23.25925 0/90 900 20.468 ±0.011 

23.26894 45/135 600 20.484 ±0.014 

23.27688 0/90 600 20.499 ±0.016 

23.28154 0/90 30 

23.28289 45/135 30 



The photometric zero point was determined assuming V = 14.582 for star A (cf. Henden et al. 2002). 
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Fig. 1. — a) An excerpt of the combined image containing tiie 0°/90° orientations. The image measures 60 X 85 arcsec. North is up and 
East is to the left. The optical afterglow of GRB 010222 (OA), Henden's bright star (A), a faint star (B) and a galaxy (G) are indicated. 
The orientations corresponding to 0° and 90° are separated by 15 arcsec in the East- West direction, b) This figure shows Q vs. U for the OA 
and G. Aperture radii of 5,6 and 7 pixels are plotted with photon noise error bars shown for 6 pixels. 



